Parametric excitation of plasma waves by gravitational radiation 
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We consider the parametric excitation of a Langmuir wave and an electromagnetic wave by 
gravitational radiation, in a thin plasma on a Minkowski background. We calculate the coupling 
coefficients starting from a kinetic description, and the growth rate of the instability is found. The 
Manley-Rowe relations are fulfilled only in the limit of a cold plasma. As a consequence, it is 
generally difficult to view the process quantum mechanically, i.e. as the decay of a graviton into a 
photon and a plasmon. Finally we discuss the relevance of our investigation to realistic physical 
situations. 
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The state of matter in regions where general relativistic treatments are desirable is often the plasma state. Neverthe- 
less, plasma physics and general relativity are quite distinct areas of physics, and accordingly there are comparatively 
few papers using a general relativistic framework that include the electromagnetic forces in their treatments. However, 
there are a few exceptions to this rule, see for example Refs. 1-5, where the plasma dynamics in a strong gravitational 
field is considered @J^], relativistic transport equations for a plasma are derived ||], a general relativistic version 
of the Kelvin-Helmholtz theorem is derived [|| and photon acceleration by gravitational radiation is considered |B. 
Omission of the electromagnetic effects for a plasma subject to gravitational forces is possible because gravity alone do 
not separate the charges, and in many cases the plasma can be treated as a neutral fluid, in spite of its electromagnetic 
properties. 

In this letter we will consider a simple model problem, that has two interesting properties: Firstly, the process of 
investigation requires a general relativistic description of the plasma dynamics to occur, and secondly we demonstrate 
the possibility of charge separation induced by the gravitational effects. We start from a monochromatic gravitational 
wave propagating through a thin plasma superimposed on a flat background metric, and consider the parametric 
excitation of a plasma wave and an electromagnetic wave. The plasma wave - that has an increasing amplitude due 
to the above mentioned mechanism - undergoes charge density oscillations. Thus the parametric excitation process 
- which of course has no Newtonian analog - shows that although no direct charge separation can be caused by 
gravitational forces, indirectly charge density perturbations may grow due to the gravitational and electromagnetic 
Q-f interaction. Theoretical aspects of our problem, i.e. satisfaction of the Manley-Rowe relations and the possibility of 
a quantum interpretation of the interaction, and the relevance of our model problem to realistic physical situations 
?-) | will be discussed in the end of this letter. 

We will refer to tensorial objects in terms of their components in a coordinate basis. Note that as observers the 
natural way to measure quantities is to project them onto a Lorentz frame {ea,a = 0, ...,3}, i.e., a frame in which 
. the metric components become those of the Minkowski metric. The frame components of the electric and magnetic 
fields E a and B a are defined by F° a = E a , F ab = S £ ^e ahc B^, where e abc are the components of the totally skew 
tensor. Expanding E = E a ea and B = B a ea to first order in the metric's coordinate components, we get that the 
electric field components in the coordinate basis are the same as in the Lorentz frame. Furthermore, the magnetic 
field components deviates only to first order in the metric components from the Lorentz tetrad components, and it is 
therefore meaningful to use the coordinate components in place of the tetrad components. This would of course not 
be the case in a strong gravitational field such as, say, the Kerr geometry. 

The Vlasov equation for the distribution / = f(x^,p a ) (where fi, v, ... = t, x, y, z and a, b, ... = x, y, z) reads Q 

where Q a = T^p^p" /p* , and this is coupled to Maxwell's equations 

1/2 



F^. v = f, f=J2^q J f(x»,p a )p»^d 3 p 



(2a) 
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*W] = - (2b) 

where p.s. stands for particle species. Here we have introduced the invariant measure Qg] 1 ^ 2 /\pt\)d 3 p on the surface 
in momentum space where g^uP^p" — —m 2 , m being the mass of the particle in question, and we use g to denote 
the determinant of the metric. The gravitational field, represented by the Christoffel symbols T^ a in the Vlasov 
equation, is assumed to be generated by some outside source. The plasma itself is assumed to generate a much weaker 
gravitational field. 

Next we assume the presence of a small amplitude gravitational wave, and we use the transverse traceless gauge. 
The line element then takes the form 

ds 2 = ~dt 2 + (1 + h+) dx 2 + (1 - h+) dy 2 + 2h x dx dy + dz 2 , (3) 

where h + = h+(u), h x = h x (u), with u=z — t, and \h+\, \h x \ <C 1. The mass of a particle with momentum p M is 
defined as m=\g liu p ti p l '\ 1 ^ 2 . From this, p l may be expressed in terms of m and the 3-momentum p a : 

p* = {m 2 + 6 abP a p b + h+ [(p x f - (p»Y) + 2h x p x py} 1/2 . (4) 
To first order in h + , h x , the Christoffel symbols read 

— r* = r* = -r x = r x =r y = — r y = -r z =r z = Ift , fKal 

xx yy tx xz ty yz xx yy 2 + ' \ UC V 

— r* — — V x — T x — —Y y —T y — —V z — !/? 

*- x y *- ty ^ yz *- tx 1 iz _ x xy 2 x ' \ ou J 

where h=dh/du. Using the Christoffel symbols (|J) of the gravitational wave, the Vlasov equation ([j]) for the unper- 
turbed (F^ — 0) distribution / = f(t,x a ,p b ) becomes 

dt p l dx a d P a [ ' 



where 



c> J = ( 1 + fr) (h+p x + h xP y 

c7"= i \ ■ j) (h x p x -h+p y 



1 

"2p* 



g*=—{ 1ABP A p B h + -2p x p y h > 



and (iab) = diag(— 1,1), A, B = x,y. Expanding p* according to p* » p* ^ — J 7 , where T = 
(h + ~/ABP A p B — 2h x p x p v ) /(2pL,), and p' ^ = (m 2 + 5 a bP a p b ) 1 ^ 2 , the Vlasov equation becomes 

dt + \ + p\ 0) ) p\ 0) dx« 9 dp- U ' [<) 

For simplicity we let the gravitational wave be monochromatic: 

h = heap [i(koz — u/ot)] + c.c. , (8) 

with the dispersion relation fco = u>o- We divide / according to / = /s,j(p a ) + fg(t, z,p a ) where the equilibrium 
distribution /sj is taken to be the Synge-Jiittner distribution and / g = / g exp [iffco 2 — wot)] ls the perturbation 
induced by the gravitational wave - which is assumed to fulfill / g <C /sj- From Eq. (0) we find 

Js — / i t \ ' 

i (w - /c p 2 /-P( ) J 

to first order in the amplitude. 

Next we assume the presence of electromagnetic perturbations with frequency and wavenumber (wi,ki) as well as 
electrostatic perturbations (w 2 ,k 2 ). The frequencies and wavenumbers are taken to satisfy the resonance conditions 



2 



luq = lui + ll>2 and ko= ki+k2. The vectors ki and k2 span a plane, and we chose the y-axis to be perpendicular to 
this plane, i.e. ki = kfx. +kfz and k2 = fcf x +&f z. Including the electromagnetic field, expanding Eq. (j?|) to first 
order in h and writing all terms proportional to h on the right hand side Eq. (m) becomes 



df n 



P a df ng 



dt dx a 



-F' 



L S bc F ac p b ' 9fns 



P\o) 



dp a 



(p'o)) 2 9x a dp" 



1 



P(o) 



[-h +lBC F aC P B + h x (F ax p y + F ay p x )} 



s.i 



dp a 



(10) 



where we have introduced / ng defined by / ng = / — / g . Furthermore, we need two components of Faraday's law. 
Expanding Eq. (2b) to first order in the amplitude h, we obtain 



d t F xy - d y F tx + d x F ty = h x (-d x F tx + d y F ty ) + h+ (d y F tx + d x F ty ) 
d t F yz - d z F ty + d y F tz = h x (d z F tx - d t F xz ) + h+ (d z F ty + d t F yz 

+h x (F tx + F xz ) - h+ (F ty + F yz ) 



(11a) 
(lib) 



In what follows, we will put h + — 0, i.e., we choose a specific polarization of the gravitational wave. It turns 
out that such a gravitational wave does not couple to a linearly polarized electromagnetic wave with magnetic field 
in the y-direction. Thus we assume the electromagnetic wave to have the opposite polarization, i.e. F yt , F xz and 
F yz are the only electromagnetic components to be different from zero in the linear approximation. Similarly F zt 
and F xt are the only nonzero electrostatic components of the field tensor. We then divide / ng according to / ng = 
fsj(.P a ) + fcm(t,p a ) exp [i(kf x + k\z - U\b)] + f cs (t,p a ) exp [i(k x x + k%z - cj 2 £)] where the time dependence of the 
amplitudes, which is due to the parametric interaction with the gravitational wave, is assumed to be slow, such that 
d/dt <C uj. In the linear approximation (i.e. no gravitational coupling) the slow time dependence vanishes, and Eq. 
© gives 



1 _ 

em 



q F yt dfi 



S.J 



iSi dp y 

qF at dfsj 
1Q2 dp a 



-ki-p/p* 



(0) 



(12a) 
(12b) 

, where the scalar 



for the electromagnetic and electrostatic perturbations respectively. Here we have uJi=Ui 

product is defined as A ■ 'B=5 a bA a B b . Considering the part of Eq. i ffoj ) varying as exp [i(fcf x + k^z — u>2t)] using 
line ar ap proximations for the factors on the right hand side we obtain / cs = + f^. The linear contribution is given 
by (12b), by replacing LU2 by W2op=^2 + id/dt (since d/dt <C u>2, division by 0J2 Op can be calculated by a first order 
Taylor expansion) and the nonlinear contribution is found to be /"J = C cs h x F yt where 



q J p x p y k 1 ■ p <9/sj ( 1 

L^os — T^T" \ —T~* 7^— + ^O^op 



1 (P( )) 3 ^i dp y 




cDi dp y 



1 /p.k a 



d 



P\o) V w i ) J d P y w iP( ) V 1 dp* 1 d P z 



(Sop/sj) 



k z lP x df S j 
P( )^i 9p z 



and C/ op = (l — p z /p\ Q \)(p v d/dp x +p x d/dp y ) + (p x p y /p^)d/dp z . Combining Eqs. (|13|) and (|2a|) we then find 



£L(u2op,k 2 )F e 



Voqh x F yt 
ik 2 



£ Jc cs d 3 P . 



where the longitudinal dielectric permittivity £l is given by ||] 



k 2 ^ J ^2op dp 



1 df SJ 3 
d p, 



(13) 



(14) 



(15) 
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(F zt ) 2 ] 1/2 - In deriving Eqs. (0) and flT|) 



p = [(p*) 2 + (p y ) 2 + (p 2 ) 2 } 1 / 2 , k 2 = [(jfcf ) 2 + (fc|) 2 ] 1/2 , and F e ^= [(i^ 4 ) 2 
we have used that the electrostatic mode is longitudinal, i.e. F xt /F zt — fcf/fcf. 

Next we turn to the part of ( fiol ) varying as exp [i(kfx + kfz — u>it)]. Again using linear approximations for the 
right hand side we obtain / cm = f\ z 



fem where fl m is given by (12a) (naturally by replacing u>i with wiop), and 



f*£. is found to be fS^ — C cm h x F es where 



r q r / fc 2 

L-em — \ kW, 



ks df, 



S.7 



lUJl 



Q 2 dp x £>2 dp z 



2 "r\ ~l~ kn — 



^0 

UJ ' 



Gopfsj 



S2 <9p x S 2 <9p z 



(16) 



Eliminating linear terms proportional to F yz and F xy with the help of Eqs. (lib) an ^ (Ha), Eq. ( J2 a| ) together with 
© gives 



^Kopjkil)^* 




where the transverse dielectric permittivity £t is given by 



£T = 1 



AW \ - 



p.s. 



LUi Op y 



(17) 



(18) 



The operators exO^op, k 2 ) and Z?(wi op , |ki|) can be expanded according to £l(w 2op , k 2 ) = £l(<^2, k 2 ) + (deL/duj 2 )id/dt 
and Z?(wi p, |ki|) = D(u>i, |ki|) + (dD/duii)id/dt. Assuming that the dispersions relations for the Langmuir wave and 
the electromagnetic wave are exactly fulfilled, we let £l(^2,&2) = D(u>i, |ki) = 0. Equations ( |l4"| ) and jl7| ) are then 
written 



QF yt 



Ci 



= C7 2 



h x F yt 
(de L /duj 2 ) 

h x F es * 
(dD/doj!) 



(19a) 
(19b) 



where the (constant) coupling coefficients Ci and C2 are trivially found from ( |14|) and ( |17[) . As can be seen, generally 
the expressions for the coupling coefficients are very complicated. In order to get transparent formulas, we present 
the result for a cold electron plasma with immobile ions constituting a neutralizing background. Letting /sj — > 5 3 (p) 
t he c oupling coefficients become C\ = C 2 = C = —ik 2x ujo/k 2 uji. Using these expressions and combining ( |l9a| ) and 
( |l9b| ) the growth rate 7 for the Langmuir wave and the electromagnetic wave is given by 



k 2x up 
k 2 ui 



U)iUJ 2 



(20) 



Since we have not included dissipation of the wave modes, the threshold value for the instability is zero. However, it is 
straightforward [jlo| to calculate the threshold by including approp riate damp ing m echanisms for the decay products. 

The fact that the same coupling coefficient C appears in Eqs. (19a) and (19b) for a cold plasma means that the 
Manley-Rowe relations Jl(J are fulfilled for that case. These relations usually follows from an underlying Hamiltonian 
structure of the governing equations, and assures that each of the decay products takes energy from the pump wave 
in direct proportion to their respective frequencies. Thus fulfillment of the Manley-Rowe relation means that the 
parametric process can be interpreted quantum mechanically - i.e. we can think of the interaction as the decay of a 
graviton with energy Hujo into a photon with energy huJi and a plasmon with energy Tiuj 2 . However, we stress tha t 
this interpretation is not always possible. When the rmal effects are taken into account, the coefficient C\ in (19a) 
is generally different from the coefficient C 2 in (|l9b|). Thus we obtain the rather surprising result that the simple 
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graviton interpretation (i.e. that the gravitational wave is built up of quanta with energy Tiujq ) of the gravitational 
field is not always applicable. 

To our knowledge, our set of governing Eqs. (jl])-(||), have not been used to consider specific physical processes 
previously. Probably this is because the system has a very high complexity - it is nonlinear and has seven independent 
variables - but, on the other hand, simpler models for the plasma cannot deal with thermal effects properly. As a 
starting point we have addressed a simple "model problem", to reduce the complexity of our governing equations. 
The main conclusion that can be drawn from our calculations are two: Firstly, wave-wave interactions can lead to 
energy transfer from gravitational to electromagnetic degrees of freedom, and vice versa. Secondly, charge separation 
and corresponding electrostatic fields may occur as a result of such interactions. However, the idealizations made 
in our calculations make it difficult to make specific physical predictions. For example, binary systems can produce 
monochromatic gravitational radiation, but the amplitude is probably to low to overcome the threshold value caused 
by plasma inhomogeneities Jllj] , except possibly in the very last collapse phase when the system produces an incoherent 
broadband gravitational pulse. Similarly, supernova events may produce high enough gravitational wave amplitudes 
for significant energy transfer to electromagnetic degrees of freedom to occur, but in most cases we expect the 
gravitational radiation to be very incoherent and broadband (there are calculations predicting highly monochromatic 
emissions of gravitational waves for certain types of neutron star formation jjjj]). Some estimates of the effects due to 
scattering of gravitational waves in supernova has been presented by Bingham et al. ||l3|| . One may also consider the 
formation of black holes as a source of gravitational radiation, especially the non-axisymmetric collapse, but - again - 
this process gives a pulse rather than a monochromatic wave. (See the review by Thorne 10 for characteristic values 
of the amplitudes, frequencies, and duration times for the processes discussed above.) 

The parametric process considered in this letter can be thought of as a gravitational analog of Raman scattering, 
where the electromagnetic pump wave has been replaced by a gravitational wave. Similarly, we can imagine that a 
gravitational analog of Compton scattering may occur. Such a process requires a kinetic treatment, and could be 
described within our formalism without very much extra difficulty. Gravitational waves may also be subject to four 
wave processes such as modulational instabilities. Furthermore, a gravitational radiation background similar to the 
microwave background is assumed to exist as a relic from big bang Jl^]. In the early universe the intensity of such 
radiation could be very high, and - if the magnitude of the gravitational waves are sufficient - processes such as those 
mentioned above could be of importance for the structure formation in the universe. It should be stated though that 
estimates of the amplitude of such primordial gravitational waves are very difficult to make, much depending on the 
insufficient data concerning the early universe. Also, much of the simplifying assumption made in this letter must be 
relaxed before such issues can be studied seriously. 
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